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Abstract 

We study the exponential Ornstein-Uhlenbeck stochastic volatility model and observe that the 
model shows a multiscale behavior in the volatility autocorrelation. It also exhibits a leverage 
correlation and a probability profile for the stationary volatility which are consistent with market 
observations. All these features make the model quite appealing since it appears to be more 
complete than other stochastic volatility models also based on a two-dimensional diffusion. We 
finally present an approximate solution for the return probability density designed to capture the 
kurtosis and skewness effects. 
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I. INTRODUCTION 



During decades the geometric Brownian motion has been widely accepted as one of the 
most universal models for speculative markets whereas nowadays we know that many ob- 
servations of real data are in clear disagreement with this model [if. A possible way out to 
these inconsistencies is to assume that the volatility is not a constant or even a deterministic 
function of time or price but a time-depending random variable. Within this approach, there 



exist the ARCH-GARCH models 



for discrete times and the stochastic volatil- 



ity (SV) models for continuous times. At late eighties different SV models were presented 
and all of them proposed a two-dimensional process involving two variables: the asset price 



and the volatility 



These and subsequent works 



lOl 111 | were basically aimed 



to study option pricing. Recently SV models have also been suggested as good candidates 



to account for ot 
dices 



ier stylized facts exhibited by the empirical observations in stocks or in- 



161 ] . Among these facts we want to single out two relevant correlations in 



time: the volatility autocorrelation and the return-volatility asymmetric correlation. 

It is well established that volatility fluctuations have long memory in the sense that the 
autocorrelation function decays very slowly with time 0, Q|. An image of daily data 
observations is the following: there is a rapid decay governed by a very short time scale (few 
days) and a fat tail with a time scale at least of the order of several hundred days. However 



common SV models result in a single exponential decay implying only one time scale 



flfl 

llJ]. Therefore we would need to add a second time scale, or even more, to properly reproduce 
the dynamics of the volatility An alternative approach to model this complex behavior 



is that of multifractal random walks which assume a power law autocorrelation 
that implies an infinite number of time scales in the evolution of the volatility 



20 



unction 



2l|. 



On the other hand, the negative return-volatility correlation -also called asymmetric 
volati lity correlation or leverage- presents a much shorter time scale of the order of 10-20 
days J22, Q, I2J]. Again, the SV models can easily incorporate the leverage effect with the 
introduction of an asymmetry parameter. This parameter also induces a bias in the return 
distribution which in turn generates the skew observed in many option smiles. Unfortunately, 
in some SV models the resulting time scale for the leverage turns out to be comparable to 



that of the volatility autocorrelation 



13. 



Q. 



Therefore, many SV models seem not to be 



able to explain the complex dynamics of the volatility. Some attempts have been made 
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to harmonize within a single model the different time scales involved 25, 26, 271. Even 
multifractal models have been recently extended to include leverage effects [2j| although the 
resulting models are quite involved and analytical results are scarce and difficult to obtain. 

We have recently addressed this problem and have extended the simple Ornstein- 
Uhlenbeck (OUl stochastic volatility model by assuming that the mean reverting level is 
itself random 25[ . We have thus considered a three-dimensional diffusion process instead of 
the usual two-dimensional one. With this model we have been able to obtain two different 
time scales, one for the leverage correlation -which, to some extend, coincides with the short 
time scale of the volatility autocorrelation- and a second scale which accounts for the long 
memory of the volatility process. The model has produced quite satisfactory analytic results 
despite the fact of the inherent difficulty in dealing with a three-dimensional process and also 
that the plain OU process has small regions where volatilit y is not positive definite. More- 
over, the OU model, either in its two-dimensional version |l5| as in the three-dimensional 
one [25), yields a Gaussian probability density function (pdf) for the stationary distribution 
of the volatility while several emp irical studies have observed that the stationary volatility 
is far from being Gaussian jll bfll Isfll]. 

All these observations have led us to consider a different SV model. Some authors de- 
fend that a good approximation for the volatility distribution is given by the log-normal 
distribution or even by the inverse Gamma distribution One of our main motivations 
here is to focus on a model whose volatility is log-normally distributed. We will thus study 
the correlated exponential Ornstein-Uhlenbeck (cexpOU) model insisting on the point of 
view of its statistical properties rather than in option pricing as have been done in previous 
studies of the model [ll|. In this way we will show how the cexpOU model, despite being a 
two-dimensional diffusion, presents a long memory behavior for the volatility with at least 
two time scales, the shortest one coinciding with that of leverage. One of the main assets 
of the model is that it is still relatively easy to handle and many analytical results can be 
obtained. 

The paper is divided into six sections. In Section [H] we present the cexpOU model 
and show its main statistical properties. Section IIHI studies the leverage and volatility 
correlations while Section IIVI applies the model to actual data thus providing a way of 
estimating the parameters of the model. The Section IS studies the probability distribution 
of the return and conclusions are drawn in Section IVII 
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II. ANALYSIS 



As mentioned in the previous section the correlated exponential Ornstein-Uhlenbeck 
stochastic volatility model has been recently addressed by Fouque et al in the con- 
text of option pricing. Here we present a slight modification of the model which consists in 
a two-dimensional diffusion process given by the following pair of Ito stochastic differential 
equations (SDEs): 

dX{t) = me Y{t) dW 1 {t) 

dY(t) = -aY(t)dt + kdW 2 (t), (1) 

where X(t) is the undrifted log-price or zero-mean return defined as 

dS IdS\ 

where S(t) is a financial price or the value of an index. The parameters a, m, and k 
appearing in Eq. ((TJ are positive and nonrandom quantities and dWi(t) = £,i(t)dt (i = 1,2) 
are correlated Wiener processes, i.e., £j(t) are zero- mean Gaussian white noise processes 
with cross correlations given by 

m)ut')) = p^(t-t'), (2) 

where pa = 1, = p {i ^ j, — 1 < p < 1). In terms of the proces Y{t) the volatility is given 
by 

a(t) = me y(i) . (3) 

For the rest of the paper we will assume that at time t the process Y(t) (and hence 
the volatility) has reached the stationary state. From Eq. |Q we see that the stationary 
expression for Y(t) is 

Y(t) = k f e- a{t - t ' ) dW 2 (t'). (4) 

This expression along with Eq. (JHJ) allow us to get the main features of the stationary 
volatility. We summarize these features as follows. 

The average value, second moment and autocorrelation of the stationary a(t) are 

(a) = me?' 2 , (a 2 ) = m 2 e 2 ^ (5) 
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FIG. 1: Representation of the daily volatility of the DJIA from 1900 to 2004 with a fit using the 
log- normal distribution given by Eq. (|9*Jl. We assume a to be the absolute value of return and a 
biassed discrepancy is clearly observed. The inset shows the same representation in a log-log scale 
in which the discrepancy is enhanced. 



and 



where 



(a(t)a(t + r)> = m 2 exp + e~ ar )} 



2c\ 



is the stationary variance of Y(t), i.e., (Y 2 (t)) = f3 2 . 

The conditional probability density function (pdf) for the volatility is 

[ln(a/m) - e'^-^ ln(a /m)] 2 



p(a,t\a ,t ) 



1 



exp 



2/3(1 - e-M*-*")) 



(6) 



(7) 



(8) 



oy27r/3(l - e - 2 "(*-*o)) 

which can be derived taking into account that a(t) is an exponential Ornstein-Uhlenbeck 
process. The stationary probability density thus reads 



(9) 
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FIG. 2: Representation of the daily volatility of the DJIA from 1900 to 2004 with a fit using the 
log- normal distribution given by Eq. l|9*|). The empirical volatility is obtained using Eq. which 
filters the biass and better fits the log-normal distribution. The inset shows the same representation 
in a log-log scale. 

We finish this section by confronting the volatility model statistics with daily data of the 
Dow Jones Industrial Average (DJIA) from 1900 to 2004. The study is limited to sustain 
the cexpOU model as a good candidate for describing the stationary volatility and we refer 



29 



3CJ). We 



tee reader to the .i_ for a m ore COmp >ete anat ysis (s ee for .stance Q[ 
first observe that from actual data we only know the time series of prices and the value of 
cr is not directly observed. In practice, one usually takes as an approximate measure of the 
(instantaneous) volatility the quantity 

\X(t + At)-X(t)\ 



a(t) 



(10) 



where X(t) is the zero- mean return. Indeed, since dX = adW then as At — > 0: 

\X(t + At) - X(t) \ « a(t)\AW(t)\ 

and this yields Eq. (ITU)) after assuming that AW(t) 2 — > At as At — > 0, but this is only true in 
mean square sense |31(. Let us check this with empirical data. Figure ^ shows the empirical 
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FIG. 3: Representation of the empirical histogram for the absolute value of returns estimating the 
volatility of the DJIA using different data sets. Circles represent the complete history of 104 years 
while diamonds represent the last 10 years. The reasonable agreement between the two data sets 
supports the assumption of stationarity. The inset shows the same representation in a log-log scale. 

histogram of the stationary volatility in annualized units j^] using Eq. (fTU|) and we see there 
a bias in real data which seems to be inconsistent with the log-normal distribution. Hence, 
the above assumptions may not be correct. An alternative breakout is to deconvolute the 
product between a(t) and the absolute value of the Wiener process 33] (see j3| for a more 
sophisticated method). Thus from the empirical data we obtain the daily return time series 
AX = X(t + At) —X(t) and, jointly with a simulation of |AW(t)|, we can derive a different 
estimation of the volatility from: 

\X(t + At)-X(t)\ 



a(t) 



(11) 



\AW(t)\ 

The result is noiser but the biass disappears as we clearly see in Fig. El We thus get a more 
reasonable agreement with the theoretical log-normal distribution assumed by our model, a 
fact that is also confirmed by more complete studies 

al. [3| have argued that the log-normal is unable to reproduce the most extreme volatility 
values and propose a distribution with power law tails. Nevertheless, this regime involves 



291 ]. On the other hand, Micciche et 
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values of sigma greater than 200% and smaller 1% (in annualized units) which may be more 
relevant for high frequency data than for the present case of daily data. 

In Fig. |2J we also plot the empirical histogram of the stationary volatility (with volatility 
estimated following Eq. ()1U|) ) taking two different data sets, one embracing the whole period 
of 104 years and a second one taking the last 10 years of data. The almost complete 
agreement of the two histograms seems to confirm the assumption of stationarity which 
is implicit in all the analysis. In any case a more positive conclusion on the stationary 
assumption would require further tests and analysis jl|, |29( . 

III. THE MAIN CORRELATIONS 

We will now present a deeper look, within the cexpOU model, on the two crucial corre- 
lation mentioned above: the leverage and the autocorrelation of the volatility. We start by 
the asymmetric volatility correlation or leverage. 

A. Leverage 

The leverage correlation is difined by 

= mX[ t + rf) t (12) 
V ; (dX{t) 2 ) 2 V ' 

and it is a measure of the correlation between the variations of return and volatility. Al- 
though the anticorrelation between volatility and return has been known for long 22 . , 



it has not been quantitatively studied until very recently when Bouchaud et al. 



24j, using 



a large ammount of daily returns for both share prices and market indices, showed that 

f -Ae- bT , if r > 0; 
£(r)= (13) 
I 0, if t < 0; 

(A, b > 0). That is, there is an exponentially decaying anticorrelation between future volatil- 
ity and past returns changes. No correlation is found between past volatility and future price 
changes thus providing a sort of causality to the leverage effect. 

We will now reproduce Eq. (fT3^) using the cexpOU model described above. From Eq. ((*[)) 
and taking into account Eq. (j3J) we have 

(dX(t)dX(t + rf) = (a(t)dWt(t)a(t + rfdW l {t + r) 2 ). (14) 



it follows from the Ito convention that if r > then dW\{t-\-r) is uncorrected with the rest 
of stochastic variables. Thus, taking into account that (dW\{t + r) 2 ) = dt, we have 

(dX(t)dX(t + r) 2 ) = (a(t)dWi (t)a(t + rf)dt (r > 0). 

On the other hand when r < the Wiener process dW\ if) is uncorrelated with the remaining 
variables. Hence, taking into account that (dWi(t)) = 0, we obtain 

(dX(t)dX(t + r) 2 } = (t<0). 

Therefore, 

(dX(t)dX(t + r) 2 ) = (a(t)a(t + T) 2 dW 1 (t))H(r)dt, (15) 
where H(t) is the Heaviside step function. In Refs. we have shown that 

(a(t)a(t + rfdWxit)) = 2pke- aT (a(t)a(t + rf)dt. (16) 

From Eqs. ©-© we get 

(cr(t)a(t + r) 2 ) = m 3 exp{2(3(e~ aT + 5/4)}, (17) 

and 

(a(t) 2 ) = m 2 e 2/3 . (18) 
Finally, substituting Eqs. (fT5]l - (fTK]l into Eq. (JT2J) yields 

C(t) = (2pk/m) exp {-ar + 2(3(e~ aT - 3/4)} H(r), (19) 

which can have the form given by Eq. (fH?|) although what one has to identify with A and b 
will depend on the different scales involved. We will discuss this issue in the next section. 
In Fig. H]we plot the leverage correlation (fTTJj) along with empirical data provided by the 
DJIA. 

B. Volatility autocorrelation 

For the cexpOU model the autocorrelation of the volatility has been already given by 
Eq. (JHJ). However, as mentioned at the end of Sec. |Hl the values of a(t) are not accessible 
from actual financial data and we have to use the instantaneous volatility in one of its 
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FIG. 4: The solid line represents the leverage correlation as given by Eq. (|19|) for the DJIA from 
1900 to 2004. Dots represent the empirical data. 

possible versions. For this reason, as a measure of the empirical volatility autocorrelation, 
we take the function C(r) = Corr [dX(t) 2 , dX(t + r) 2 ], that is, 

r( , = (dX(t) 2 dX(t + r) 2 }-(dX(t) 2 } 2 
1 ' {dX{t) A ) - (dX(t) 2 } 2 ' 1 J 

Again in Ref. j^j], we have shown that this correlation function can be written as 

C{T) = 3<<r*> - (a(tW ' (21) 

The averages appearing on the right hand side of this equation can be evaluated using 
Eqs. flBJ-Q- The final result reads (see Fig. |3J) 

exp[4/3e" ar ] - 1 



C(r) 



3e^ - 1 



(22) 



IV. DETERMINATION OF PARAMETERS AND MULTIPLE TIME SCALES 

In the cexpOU model, there are four parameters to be determined: a, j3 (or k), p and 
m. We estimate these parameters as follows. Firstly, from the long time behavior of the 
volatility autocorrelation (see Eq. and Fig. El below) we can determine the value of 

10 



a. Furthermore, from the short time behavior of C(r) we can estimate k and hence (3 (cf. 
Eq. (|27jl anf Fig. [Hjbelow). After knowing (3 the value of the normal level m can be obtained 
through the return variance, that is (see Eqs. © and (fTU)l). 

(AX 2 ) « mV^At, 

where At = 1 day. Finally, the correlation coefficient p is estimated from the value of the 
leverage correlation, £(r), as r — > + , that is (see Eq. (fTU|) and Fig. HJ), 

£ ( o+) = ^V/ 2 . 

m 

For the DJIA and using daily data from January 1, 1900 to June 14, 2004 (a total of 
28,540 data points) the estimated parameters are 

a = (1.82 ±0.03) x 10~ 3 days -1 , k 2 = (1.4 ± 0.2) x 10~ 2 days" 1 , (23) 

which implies (3 = 3.8 ± 0.6. Finally 

m = (1.5 ± 0, 4) x 10" 3 days~ 1/2 and p = -0.4. (24) 

Note that the value of m corresponds to an annual normal level of volatility of (2.4 ±0.7)%. 



A. Time scales 

Once determined the relevant parameters we can address the crucial question of how 
many time scales have our model. Recall that the empirical volatility needs, at least, two 



time scales for dail y d ata |l£ 
models 

Eq. ((22) in tli e form 



In contrast with other two-dimensional stochastic volatility 
the cexpOU model has this feature. Indeed, we can write 

C(t) = —r^ Y ^E2- e - na \ (25) 

which indicates that there are infinite time scales. In Fig. Owe plot In C(r), with C(r) 
given by Eq. (|25j) . as a function of ar for different values of the parameter (3 = k 2 /2a. We 
there observe two relevant time scales as limiting cases of an infinite set of scales (see below). 
Moreover, the bigger (3 is, the more distant the two time scales are (see Eq. (J2*%j) and Fig. |3J). 
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FIG. 5: We plot in C(r) given by Eq. (|25|) as a function of ar for different values of j3 = k 2 /2a 
in a semi log-scale. The figure clearly shows the existence of two asymptotic time scales separated 
by a sum of multiple time scales. This effect is enhanced in bigger values of (5. 

Had we we plotted other two-dimensional SV models we would have 

a straight line unable to fit the curve observed in real data. 

Let us go deeper in this issue. Suppose first that we are looking at C(r) at long times 
such that ar ^> 1. In this case, neglecting trascendentally small terms of the order of e~ 2ar 
and higher in Eq. ()25|) . we have 

C(r) « Jf e~ ar . (26) 

V > 3e 4/3_ 1 V ) 

Thus, the long time behavior of the volatility autocorrelation is governed by the characteristic 
time a~ l (see Fig. |HJ). On the other hand, for short times such that ar <C 1, one can see 
from Eqs. (I7j) and lf22) that 

C(r) = — l — fe 4 ' 3 - 2fcV - ll + O (a 2 r 2 ) (27) 



Ze^ - 1 

(see Fig. El). In this case the short time behavior of Cir) has the characteristic time l/2k 2 . 
In many situations a -1 is of the order of few hundred days (for the DJIA = 549 ± 9 
days) while l/2k 2 is of the order of few weeks (for the DJIA l/2k 2 = 35 ± 18 days) which is 
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FIG. 6: The solid line is the volatility autocorrelation as given by Eq. (|22[) and the dashed lines are 
the asymptotic representations of this equation. Dots represent the empirical data for the DJIA 
from 1900 to 2004. 
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FIG. 7: Autocorrelation of the process Y(t) in log-log scale. The solid line is the autocorrelation 
of y = ln(<r/m) for the DJIA from 1900 to 2004. Dots represent the empirical data. 
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consistent with the separation of time scales just mentioned. Therefore, within the cexpOU 
model, the autocorrelation of the volatility presents, at least, two time scales: the longest 
one governed by Ti ong — 1/a and a shorter one governed by r s h OTt = l/2k 2 . We observe that 
the dimensionless parameter j3 defined in Eq. (J7J) can be written in terms of Ti ong and r s h ort 
as 

P=-P^, (28) 
which shows that /3 is a measure of the distance between the long and the short time scales 
(see Fig. [5J). 

We can test the soundness of the model by looking at the log- volatility autocorrelation, 
that is, the autocorrelation of the random process Y(t) (see Eq. (J3J)). One can easily see 
from Eq. (jlj) that the cexpOU model has an exponential time decay for the log-volatility: 

Corr[F(t + r)Y(t)] = {k 2 /2m)e" aT . (29) 

Let us thus see whether the data present the same behavior. To this end we start from 
the absolute value of the return as an estimator of the daily volatility (cf. Eq. (fTU|) with 
At = 1 day) we then take the logarithm of this quantity and evaluate its autocorrelation. 
We finally compare this empirical correlation with the theoretical one given by Eq. (|2*9"|) . 
In Fig. we show both correlations in log-log scale. We observe that in plotting Eq. (|29|) 
we have used the same parameter values as in Figs. 0] and El The agreement between the 
theoretical model and the empirical data is quite remarkable except, perhaps, at very short 
time lags where one can observe a small but systematic deviation 35]. We can therefore 
conclude that the cexpOU model provides a fairly satisfactory description of the multiscale 
behavior of the volatility. 



B. The leverage 



The same kind of asymptotic analysis we have performed on the volatility autocorrelation 
can now be applied to the leverage correlation (fT9*|) . Thus for long times, ar ^> 1, we have 
(r>0) 

(30) 



CM = -pke-^ 2 \e~ aT + O (e~ 2aT ) 
m 1 v ' 

which has the form of Eq. (JT3J) with b = a and 

A = -pke-^ 2 = A 
m 



long- 



(31) 
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However, as we can realize looking at Eq. (|28|h the value of (3 will be usually large, at least 
larger than 1. This is precisely the case of the DJIA where (3 ~ 3.8 > 1. Consequently Ai ong 
is exponentially small and the long-time behavior of £(r), as expressed by Eq. (|3U|) . turns 
out to be neglegible and practically undetectable in empirical observations. Therefore, the 
effect of the leverage correlation has to be sought in the short-time regime ar < 1. In this 
situation an expansion of Eq. (fT^j) yields 

CM ~ -pkeP /2 e- k2r (or < 1), (32) 
m 

where we have taken into account that a <C k 2 . Again Eq. ()32|) has the form of Eq. ()13|) 
with b = k 2 and 

A = -pke^ 2 = A short . (33) 
m 

Note that in those empirical situations in which /3 ^> 1: 

-**-short 

thus confirming that the long-time behavior of the leverage correlation is neglegible and 
practically undetectable. Therefore, the leverage correlation is governed by Eq. (J32|) which 
implies one single time scale given by r/ e „ = 1/k 2 . Note that this time scale is of the same 
order of magnitude (a litle longer) than the short time scale of the autocorrelation of the 
volatility (cf. Eq. flZZj)). 



V. THE DISTRIBUTION OF RETURNS 



We will now address the question of the distribution of returns chracterized by the return 
pdf, p(x, t), or, in an alternative and equivalent way, by its characteristic function ip x (uj,t). 

Let P2{x, y, t\yo) be the joint probability density function of the two-dimensional diffusion 
process (X(t),Y(t)) described by the pair of stochastic differential equations as given in 
Eq. (0). This pdf obeys the Fokker-Planck equation 

dp2 d l l2 d 2 P2 .j d 2 1 2 2 d 2 p 2 ( v 

■dt = a dy- (yP2) + 2 k W + Pkm d^y- (e P2) + 2 m * ~W> {M) 

with initial condition 

P2(x,y,0\yo) =S(x)5(y-y ). (35) 
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It does not seem to be possible obtaining an exact solution to this problem and therefore 
we will search for approximate expressions. There are several possible strategies to this end 
depending on the particular values of the parameters of the model for a given market. Here 
we will treat the case in which the "vol of vol" k is much greater than the normal level of 
volatility exemplified by m. In other words we will assume that the parameter 

A = - » 1 (36) 
m 

is large. Note that this is the case of the DJIA daily data, since from Eqs. (}2*3*|) and we 
see that A ~ 10 2 . 

We define a dimensionless time, r, and two scaling variables, u and v , by 

r = k 2 t, u = Xx, v = Xy. (37) 

Then the FPE (El reads 



d P2 1 9/ v, i^d 2 ^ d 2 / v/x n 1 2v/x d 2 p 2 . . 

and 

p 2 (u,v,0\v ) = 5(u)8(v - v ). (39) 

We can remove the dependence on the initial volatility of the joint distribution by performing 
the average of p 2 {u, v,t\vq) over all possible values of vq. If we assume that the initial 
volatility is in the stationary state, then the unconditional joint density is given by 



p(u,v,r)= I p 2 (u,v,T\v )p st (v )dv 0: 

where 



1 



P*t( v o) = /0 a ,, e 



Hence, p(u,v,r) is the solution to the initial-value problem: 

dr ~ 2(3dv [VP) + 2 A dv* + Xp dudv ^ V > + 2 & du^ (4 ° } 

and 

p(u,v,0) = -^=^e-^5(u). (41) 
We can write a more convenient equation for the characteristic function defined by 

/oo rco 
e luJlU du / e iuJ2V p(u,v,T)dv. 
-oo J — oo 
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The Fourier transform turns problem (|4U( ) -(|4i p into 



dip 



and 



1 d f X \2 2 / \ 

Xpu 1 uj2ip (u)x,W2 - i/X, t) - 2^1^ ( u u u 2 - 2i/A, r) 



^(^,o; 2 ,0)=e-^/ 2 . 



(42) 



(43) 



Once we know the solution to this problem, the marginal charactistic function of the return 
(p x ( UJ i, T ) is obtained through 



(44) 



In the Appendix A we show that an approximate expression for fx (<^i > T ) > valid for large 
values of A, is given by 



fx (u/X,r)^e-^ 2x2 



3 4 

1 - Upf3 2 a(T/2{3)^ + 4/3 3 6(r/2/3)^ + 0(1/A 5 ) 



(45) 



where the functions a(z) and b(z) are defined in the Appendix A by Eqs. (|A10|) and (jAllj) 
respectively. 

The inverse Fourier transform of Eq. (|43|) finally results in an approximate expression 
for the pdf of the return which in the original variables (cf. Eq. (pTTj) ) reads 



Px{x,t) 



\Jl-KmH 



-x 2 /2m 2 t 



1 - 



pka(at) ( x 
x 



k 2 Hat) 

+ „ , — -"4 



(46) 



Mat) 2 *\ y^H, 

where H n (x) are Hermite polynomials. Observe that Eq. (|46() is an expansion that corrects 
the Gaussian density: 



p(x,t) 



1 



-x 2 /2m 2 t 



(47) 



\Zz7rrnH 

This density would correspond to the return pdf if the volatility would have been a de- 
terministic quantity, i.e., if Y(t) = 0. We incidentally note that Eq. (J4lj|) constitutes an 
example of Edgeworth series, a well stablished approximation procedure in mathematical 
statistics 2f| • 

The deviation from the Gaussian character of the return pdf as shown in Eq. (146(1 is also 
evidenced by the existence of non-zero cumulants of order higher than two, and specially 
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FIG. 8: The probability distribution of 20 day return for the DJIA. We also represent the approx- 
imate density presented in Eq. (jlHjl toghether with the Gaussian pdf given by Eq. (|37|) .The inset 
shows the same representation in a semi-log scale. 

by the skewness, 71, and the kurtosis 72. These are related to cumulants by 7! = k 3 /k^ 2 
and 72 = K4/K2. Recall that, knowing the characteristic function, cumulants can be easily 
obtained by 

Qn 



uj=0 



Finally from Eq. (}45|) we get 



f— a(at) „b(at) , , 

^W^i. («) 

In Figs. |H1 and 03 we represent the empirical probability distribution of 20 day returns 
for the DJIA. We adjust the empirical distribution with the approximate density given by 
Eq. which, as shown in these figures, fits rather well the empirical density in any case 
better than the Gaussian pdf (|37j) also represented. 
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return absolute value 



FIG. 9: The positive and negative tails of the probability distributions reported in Fig. |H1 
VI. CONCLUSIONS 

We have revisited the exponential Ornstein-Uhlenbeck stochastic volatility model. Pre- 
vious studies on the model were mainly focussed on option pricing somewhat overpassing 
the description of the underlying asset We have shown that the model explain fairly 
well some relevant properties of the market. A good feature is that the model has a log- 
normal stationary probability distribution for the volatility consistent with data. We have 
found that the model can also explain the observed long-range behavior (with at least two 
time scales) of the volatility autocorrelation providing also an infinite cascade of time scales. 
Moreover, we have seen that the values of the parameters that adjust the empirical volatility 
autocorrelation are also the appropriate ones to describe the short range return-volatility 
correlation. This consistency has been checked for the particular case of the daily Dow Jones 
index from 1900 to 2004. 

Perhaps one of the main reasons to consider the cexpOU model lies on its simplicity, 
without the need to resort to fractal phenomena. It still keeps a two-dimensional diffusion 



formalism which is fully ^e^crjbec 
models in the literature 




by only four parameters. Other two-dimensional diffusion 
10Lll5lll6j| cannot describe the multiple time scale behavior 
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and they need to add a third diffusion equation with a second time scale in order to account 



for the long memory of the volatility [25|, |27[. The only drawback of the cexpOU is that 
we do not have an exact solution for its characteristic function in contrast with other SV 
models jisil"!^. However, we have also presented an approximate solution starting from the 
standard Black-Scholes Gaussian distribution. The solution takes into account higher order 
cumulants and can easily incorporate the effects of skewness and kurtosis. The resultant 
distribution improves the Gaussian pdf in the statistics of the Dow Jones 20 day returns 
and may be interesting for option pricing. 
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APPENDIX A: THE CHARACTERISTIC FUNCTION OF THE RETURN 

In order to prove Eq. (|4"5j) we will look for an approximate expression of the joint distri- 
bution (f(uJi, ct>2, t) valid for large values of A. We first note that the marginal charateristic 
function of the return can be obtained from the joint characteristic function by setting uj2 = 
(see Eq. (|44|) ) . Therefore, we will look for a solution to the problem (j42|) - (j4*3*j) that for small 
values of U2 takes the form: 

(p(ui,U2,r) — exp j— A(<jJi,t)u% + B(ui,t)u2 + C(ui,t) + 0(u;f) J . (Al) 

Substituting this into Eq. (J4*2~j) yields 



Au 2 2 + Bu 2 + C = -^-u 2 {2Auj2 + B) + \\ 2 uj 2 2 
[A [2icj 2 /A+ 1/A 2 ] +iB/\} 
+ -u^exp [AA [iu 2 /X + 1/A 2 ] + 2iB/\} , 

where dots represent a time derivative. Expanding the exponentials in powers of 1/A up to 
first order and equating equal powers in u 2 , we obtain an approximate first-order system 
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of ordinary differential equations for A, B and C with initial conditions A(ui,0) = j3\ 2 /2, 
B(ui, 0) = and C(u\, 0) = 0. The solution to this problem is straightforward and reads 



A{w u t) ~ ^A 2 e" 2 ^ l)r + A 2 [l - e~ 2 ^ l)r ] /4 7 (o;i), 



(A2) 



B^, r) ~ i\u 2 [(2/3 - 1/7(^1)) e-^ T + (l - e -^) T ) /7 2 (^i)] (l - e~^ T ) /2 7 (^) 
+ Apw^l-e-^^^wi), (A3) 

CVx, r) ~ ^ 2 r/2 + Wl 3 [r - (l - e~^ T ) / 7 (^)] / 7 (^) 

r + (l - e ~ 2 ^ T ) /2 7 ( Wl ) - 2 (l - e^^) T ) / 7 K) 



- w 



+ /3(l- e - 7(wi)T )' 



/2 7 2 K), 



where 



(A4) 



(A5) 



From Eqs. (|44|) and (jAl|) we see that the marginal characteristic function of the return 
reads 

<px(wuT) = exp[-C(wi,T)]. (A6) 

Now the return pdf will be given, in the original variables (cf. Eq. (|37jl ). by the inverse 
Fourier transform 

Px(x, t) = ±- r e-^ipx (w/A, r) du. (A7) 

Z7T J— oo 

We write fx (^/A,r) in the form 

(<*V A, r) ~ exp {-cu 2 r/2A 2 - ipu 3 Z(w/\, r)/A 3 + wV(w/A, r)/A 4 } , 

where 

£(w/A, r) = [r - (l - e^^) / 7 ( W /A)] / 7 (w/A), 

and 



i/;(uj/X,t) 



t + 1 - e 



- 2 ^ /a)t )/2 7 (cu/A) - 2(1 - e^ /A)T )/ 7 (c^/A) 



+ ^fl_ e -2 7 (./A)r 



/2 7 2 (^/A) 



Taking into account the expression for j(u/X) given by Eq. ()A5|) and expanding £ and "0 in 
powers of 1/A, we get 

Vx ( W /A, r) ~ e- 2 ^ 2 exp |-4^ 2 a(r/2/3)^ + 4/? 3 6(r/2/5)^ + 0(1/A 5 )| (A8) 
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which after a further expansion can be written as 

Vx (uj/X,r)^e-^ 2X 
where 



3 4 

1 - MpP 2 a(r/2P)^ + 4/3 3 6(r/2/3)^ + 0(1/A 5 ) 



and 



a(z) = z — (l — e , 
b(z) = (1 + 2p 2 )z + 2p 2 [ze- z - 2 (l - e~ 2 )] - (l - e~ z 



(A9) 



(A10) 
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